Popular Summary 

Persistent Nature of Secondary Diurnal Modes of Precipitation 
over Oceanic and Continental Regimes 
Song Yang, Kwo-Sen Kuo, and Eric A, Smith 

The first joint Tropical Rainfall Measuring Mission (TRMM) by the US National Aeronautics 
and Space Administration (NASA) and the japan Aerospace Exploration Agency (JAXA) 
launched on November 22, 1997 provides unique long-term precipitation measurements over the 
global tropics. This study develops a powerful clock-face type graphic display method to dearly 
and accurately show the global patterns of the rainfall diurnal phase and amplitude at the same 
time, using the 8-vear TRMM rainfall datasets. The late evening early morning diurnal phase 
over (Kean and mid-late afternoon diurnal phase over land are graphically displayed consistent 
with the global continents. Ibis paper is also the first to show the globally spatial patterns of the 
minor diurnal phase and intensity of rainfall diurnal cycles. The impact of horizontally spatial 
scales on rainfall diurnal cycle is discussed. Results show that any selected spatial scale has no 
impact on the major diurnal peak of rainfall, but a significant impact on the minor diurnal peak. 
The advantages and disadvantages of using a harmonic analysis method in studying the rainfall 
diurnal cycle are discussed in detail. This method can exhibit global spatial patterns of the major 
diurnal peak, while suppressing the influence of rain noise on diurnal cycles. However, it reveals 
a weaker diurnal intensity and shifted phases. The possible false distribution of the minor 
diurnal peak caused by this method is physically and graphically explained. 
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Abstract 


This investigation seeks a better understanding of the assorted mechanisms controlling the 
global distribution of precipitation diurnal variability based on the use of Tropical Rainfall 
Measuring Mission (TRMM) microwave radiometer and radar data. The horizontal distributions 
of precipitation's diurnal cycle arc derived from eight years of TRMM Microwave Imager (TMI) 
and Precipitation Radar (PR) measurements involving three TRMM standard rain rate retrieval 
algorithms — the resultant distributions analyzed at various spatiotemporal scales. The results 
reveal the prominent and expected late-evening to early-morning (LE-EM) precipitation maxima 
over oceans and the counterpart prominent and expected mid- to late-afternoon (MLA) maxima 
over continents. Moreover, and not generally recognized, the results reveal a widespread 
distribution of secondary maxima occurring over both oceans and continents — maxima which 
generally mirror their counterpart regime's behavior. That is, many ocean regions exhibit clear- 
cut secondary MLA precipitation maxima while many continental regions exhibit just as evident 
secondary LE-EM maxima. This investigation is the first comprehensive study of these globally 
prevalent secondary maxima and their widespread nature, a type of study only made possible 
when the analysis procedure is applied to a high-quality global-scale precipitation dataset. 

The characteristics of the secondary maxima are mapped and described on global grids using 
an innovative clock-face format, while a current study to be published at a later date provides 
physically-based explanations of the seasonal-regional distributions of the secondary maxima. In 
addition to an “explicit" maxima identification scheme, a “Fourier decomposition - ’ maxima 
identification scheme is used to examine the amplitude and phase properties of the primary and 
secondary maxima — as well as tertiary and quaternary maxima. Accordingly, the advantages, 
ambiguities, and pitfalls resulting from use of Fourier harmonic analysis are explained. 





1. Introduction 


Studies of atmospheric diurnal processes which are influenced by the regulated daily 
cycle of incoming solar radiation at the top-of-atmosphere {TO A) have been taking place for 
over one hundred years. The seminal study of Harm (1901) was the first to address 
precipitation's diurnal cycle. Observational and modeling analyses have demonstrated that 
diurnal processes are evident in many atmospheric quantities. These include precipitation (e.g.. 
Hong et al. 2005 and Yang and Smith 2006), surface temperature (e.g.. Smith 1986), surface 
winds (e.g., Deser and Smith 1998), surface pressure (e.g.. Petenko and Argentini 2002). vertical 
motion (e.g.. krishnamurti and Kishtawal 2000). cloudiness (e.g.. Wylie and Woolf 2002), and 
surface (e.g.. Smith et al. 1986) and TOA (e.g.. Smith and Rutan 2003) radiation fluxes — as 
foremost variables. Notably, lengthy time series derived from satellite measurements of various 
atmospheric variables have motivated new types of research concerning diurnal variability. 

This has been particularly true for precipitation since the advent of Tropical Rainfall 
Measuring Mission (TRMM) and its associated TRMM Microwave Imager (TMI) radiometer 
and Precipitation Radar (PR) rain rate retrievals dating back to November 1997. Many 
mechanisms have been proposed to explain precipitation diurnal behavior. The “static radiation- 
convection” (SRC) mechanism (e.g.. Randall et al. 1991), "dynamic radiation-convection” 
(DRC) mechanism (e.g.. Gray and Jacobson 1977). “nighttime radiative cooling driven elevated 
relative humidity” (NRC-ERH) mechanism (e.g.. Tao et al. 1996; Sui et al. 1997. 1998). or the 
slowly evolving, diurnally-varying "large scale vertical motion" (LS-VM) mechanism (e.g., 
McBride and Gray 1980) represent various possible explanations for the well-established "late 
evening-early morning” (I. E-EM) oceanic surface precipitation maximum, while the diurnally- 
regulated surface solar radiative heating (SRH) mechanism, which can manifest itself in either 
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the form of a "static destabilization” (SD) mechanism or a “differential heating” (DU) 
mechanism, are the foremost explanations for the often-observed “mid- to late-afternoon” 
(ML A) continental surface precipitation maximum (e.g., Ramage 1971 and Pielke 2002). A 
number of studies using regional observations (e.g.. Schwartz and Bosart 1979; Oki and Musiake 
1994: Anderson el al 1996; Sui et al 1998; Chen et al. 1999; and Dai 2001) have reported a 
secondary LE-EM peak in the continental precipitation diurnal cycle. The main possible 
explanations for this mode, which have been reviewed by Yang and Smith (2006). consist of the 
“mobile terrain-forced precipitating system" (MTFPS) mechanism and a continentally-based 
NRC'-ERH mechanism. There is also a counterpart MLA secondary peak in the oceanic 
precipitation diurnal cycle (e.g.. McGarry and Reed 1978; Reed and Jaffe 1981; Augustine 1984: 
Fu et al. 1990: Serra and McPhaden 2004), for which the only plausible explanation that can be 
put forward is an ocean surface heating (OSH) mechanism, which can involve the near surface 
water layer (for oceans, seas, and inland lakes) and/or the moist boundary layer over the water. 

However, there is yet to be and investigation of the secondary diurnal mode of 
precipitation as a global phenomenon. Yang and Smith (2006) provide a detailed review and 
analysis of the variety of mechanisms controlling the diurnal cycle of precipitation. Based on the 
use of TRMM precipitation data, their analysis demonstrated that precipitation diurnal variability 
is a global phenomenon with embedded diurnal forcing factors — far more complex than can be 
explained with a few general causes (i.e.. the approach followed by nearly al! past literature 
concerning this topic). The primary and dominant LE-EM oceanic precipitation maximum is 
often accompanied by a secondary MLA maximum, while the primary and dominant MLA 
continental precipitation maximum is often accompanied or even replaced by a secondary LE- 
EM maximum. 
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As Yang and Smith (2006) emphasize, there are a host of mechanisms at work that 
produce the precipitation diurnal process, not one of which solely explains the entire process. 
Instead, a mixture of mechanisms whose individual components control regional and smaller 
scale diurnal modes work together to produce the averaged global process. Notably, most 
modeling studies investigating diurnal variability, particularly studies based on the atmospheric 
general circulation model (GCM). have focused on stand-alone monolithic mechanisms (e.g.. 
Randall et al. 1991, Dai and Trenberth 2004). The recent study of Dai (2006). reviewing results 
from 18 coupled atmosphere-ocean GCM (AOGCM) models, indicates large discrepancies 
involving diurnal variability of precipitation amongst the different models. Therefore, given 
these findings and our experience with analyzing TRMM data, it appears that modelers should 
improve the "diurnal processes' capabilities of their models — including precipitation — to 
ensure that they reliably simulate atmospheric circulations that themselves are diurnally 
modulated. 

Furthermore, more detailed observational analyses of the spatiotemporal characteristics 
of diurnal cycle variations are needed to reveal their unique properties so that the mechanisms of 
the diurnal cycle can be better physically explained. Modelers would then have a better 
quantitative basis from which to refine model processes that control the diurnal variability of 
precipitation. Yang and Smith (this issue) describe the seasonal climatology of the diurnal cycle 
at global scales, its seasonal variability, and the contrasting behaviors of convective and 
stratiform components — demonstrating that the secondary diurnal mode is largely modulated by 
the diurnal cycle of stratiform precipitation. In addition. Yang et al. (2006) found consistent 
appearance of global scale secondary’ precipitation maxima from eight years of TRMM data. 
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As a follow-up. this study demonstrates the widespread consistency in the occurrence of a 
secondary diurnal mode in precipitation globally, and at seasonal and regional scales. To obtain 
the detailed view, an innovative clock-face graphic display scheme is used to aid the analysis. 
Eight-year TRMM precipitation datasets derived from three level 2 TRMM standard rain rate 
algorithms are used to conduct the analysis. A Fourier harmonic analysis scheme is also applied 
to highlight the dominant diurnal modes of precipitation variability while suppressing the 
unimportant diurnal harmonics. The advantages, ambiguities, and pitfalls related to use of 
harmonic analysis are also discussed to draw attention to what can and cannot be achieved by 
using Fourier decomposition analysis in study ing precipitation's diurnal cycle. 

2. Methodology and Dataset 

An "explicit" precipitation maxima selection scheme is first used to identify the primary 
and secondary diurnal modes at various space-time scales. A "Fourier decomposition" selection 
scheme is then applied as a cross-checking method and as a means to filter out ambiguous time 
series noise, (iridded properties of the diurnal cycle are illustrated using the clock-face display 
scheme. The different colors of the clock-hands denote the different diurnal modes, while the 
length and position of the clock-hands denote amplitude and phase characteristics of the separate 
diurnal modes. In addition, clock-face-colors arc used to express the canonical phase interval of 
the primary diurnal mode. 

The existence of a secondary peak in a diurnal sequence of surface precipitation rate is 
found using two different precipitation maxima selection schemes. The first is referred to as the 
“explicit" selection scheme based on logic-based detection of local relative precipitation rate 
maxima in a given diurnal time series. This scheme consistently identifies the appearance of 
either primary or secondary modes, as long as each of the separate amplitudes exceed a 
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minimum amplitude-emergence threshold and as long as there is sufficient separation in the 
positions of two detected peaks according to a second phase-angle separation threshold. 
Application of the "explicit" scheme may result in identification of both primary and secondary 
modes, only a primary mode, or no modes at all. Tertiary and quaternary modes are also 
identified with the "explicit” scheme. 

The second scheme is referred to as the "Fourier decomposition" selection scheme. This 
scheme by definition identifies amplitudes and phases of all wave numbers. However, in the 
analysis procedure we apply, only the primary' through quaternary modes are examined. It is 
noted that assignment of phases for the secondary through quaternary modes requires an 
adjustment procedure to the raw phase information for wave numbers 1 - 3 (i.e., Fourier 
harmonics 2 - 4). For example, in assigning the phase of a secondary mode, and defining the raw 
phase angles of the 1 M and 2 nd Fourier harmonics as 0, and <j > 2 . existence of a peak for the 
secondary diurnal mode is only established when each of the following two conditions is 
satisfied: 


\<t> 2 -0,| <45° 

(1) 

|n + <t> 2 ~ 0,| < 45° 

(2) 


Thus, the peak of the 2 lid harmonic is separated from the peak of the I s ' harmonic by at least six 
hours. Otherwise, the primary and secondary modes would merge into a different primary mode 
with a different phase. Likewise, the phase positions of the tertiary and quaternary modes must 
be established by assuring appropriate phase-angle separations. As indicated earlier, the clock- 
face graphic display scheme can then be used to illustrate the various parameters associated with 
the diurnal analysis on a regular grid. 
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Eight-year (1998-2005) precipitation datasets derived from the three main level-2 TRMM 
standard precipitation algorithms are used for the study. These are the TMl-only (i.e, alg 2a 12). 
PR-only (i.e., alg 2a25), and combined PR-TMI (i.e., alg 2b31) rain rate retrieval algorithms. 
The level-2 instantaneous rain rates at native orbit-swath pixel locations are first binned into 
eight (8) 3-hourly mean solar time (MST) bins — then grouped at different spatial and temporal 
resolutions. A hierarchy of 5°x5°, 10°xl0°. 20°x30° and 20°x60° latitude-longitude spatial grid 
resolutions are then used in processing the data at seasonal and 8-year mean temporal scales. 

3. Diurnal Properties of Precipitation at Global Scale 

The distribution of the 8-year mean TRMM-based precipitation diurnal cycle over oceans 
and continents depicts a distinct climate feature on the global scale. Figure 1 shows the climate 
characteristics of diurnal variability from TMl-only, PR-only, and combined PR-TMI rain rate 
retrievals. This diagram clearly shows that the precipitation diurnal cycles from the three 
TRMM rain rate products are consistent in terms of amplitude and phase over oceans, and 
consistent in phase over continents although with somewhat larger TM1 precipitation amplitudes. 
Oceanic precipitation has a dominant late evening maximum in the 3-6 MST timeframe while 
continental precipitation has a prevailing late afternoon peak in the 15-18 MST timeframe. 

Seasonal variability of the diurnal cycle is another important feature described in detail 
by Yang and Smith (2007). Figure 2 illustrates the climatology of seasonal variations of the 
diurnal cycle over oceans and continents using 8-year 2a25 PR rain rates. The primary 
maximum at 3-6 MST of oceanic precipitation is prominent in Spring, Autumn and Winter, but 
relatively weak in Summer. At the global scale, the secondary maximum over oceans is not 
obvious for any season. The primary late afternoon peak at 15-18 MST of continental 
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precipitation is the most dominant feature. The secondary late evening maximum is most 
apparent in Summer and Winter, while weak but detectable in Spring and Autumn. 

The diurnal variations of convective and stratiform precipitation from 2a25 are also 
illustrated in Figure 2. Over oceans, convective precipitation exhibits the same diurnal 
properties as total precipitation, while stratiform precipitation exhibits a secondary peak in early 
afternoon in Spring and Winter in addition to its dominant late evening maximum. Over 
continents, the dominant late afternoon afternoon peak in convective precipitation mimics the 
diurnal behavior of total precipitation, while stratiform precipitation exhibits a late evening 
maximum at 3-6 MS I and a late afternoon maximum at 15-18 MST. The TMI-only and 
combined PR-TM1 rain rate products (figures omitted) indicate mostly similar characteristics 
insofar as diurnal variability. 

4. Persistent Nature of Secondary Diurnal Mode 
4. / Characteristics of Secondary Mode 

Horizontal distributions of the precipitation diurnal cycle based on the 8-year TRMM rain 
rate products analyzed on a 5°x5° grid are shown in Figure 3 illustrated in the format of clock- 
face grids. A blue (green) color-face of a given clock denotes a pre-noon half-day (post-noon 
half-day) phase interval of the dominant maximum for each grid location, where the pre-noon 
half-day phase interval is defined from 0-12 MST and the post-noon half-day phase interval from 
1 2-24 MST. The transition of clock-face color is evident along almost all coastlines. This draws 
attention to the dominance of the pre-noon precipitation maximum over oceans (blue) and the 
post-noon precipitation maximum over continents (green). This behavior was reported by Yang 
and Smith (2006) for a one-year dataset ( 1998). Here it is graphically and robustly illustrated 
using 8 years of data at high spatial resolution at specific grid regions. It is seen that at the 5°x5° 
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spatial scale, a few oceanic grid locations exhibit dominant post-noon maxima, while a few 
continental grid locations indicate dominant pre-noon maxima. The reasons for these exceptions 
deserve future examination. In addition, close examination of the clock-faces reveals that 
secondary maxima (denoted by white inner faces) are widespread throughout the tropics and sub- 
tropics observed by the TRMM satellite, in which secondary post-noon peaks are found over 
oceans and secondary pre-noon peaks found over continents. These characteristics are highly 
consistent between the TMI. PR. and combined PR-TMI rain rate products. 

Similar analyses conducted at 20°x30° spatial resolution scale are illustrated in Figure 4. 
At this grid scale, there is no ambiguity as to the separation of the primary' pre-noon and post- 
noon modes over oceans and continents, respectively (with the proviso that the mixed land-water 
Maritime continent exhibits purely oceanic diurnal precipitation properties in regards to the 
primary mode). Note that at this grid scale, .there is perfect agreement amongst the three 
algorithms insofar as the arrangement of the blue and green clock-faces. Also, as in Figure 3, 
there is widespread but not complete distribution of the secondary mode (marked by white inner 
clock-faces), noting that for this diurnal property, whereas there is general agreement amongst 
the three algorithms, there is not perfect agreement. These two aspects of the secondary mode, 
along with other more obscure characteristics found by close examination of the amplitudes and 
phases of the secondary mode (i.e., denoted by the appearances, lengths, and positions of green 
clock-hands) indicate that it is a more complex process than that of the primary mode. It is also 
apparent that the secondary mode over continents is more prevalent. Quantitatively, the results 
show that there is 72%. 61%, and 77% consistency in the occurrence of secondary peaks 
identified between the three algorithm pairings 2al2 - 2A25. 2a 1 2 - 2b31, and 2a25 - 2b31. 
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Figure 5 presents another set of similar analyses, but now based only on the PR algorithm 
and comparing the 20°x30° result shown in Fig. 4 with results at I0°x!0° and 20°x60° spatial 
resolutions. Although there is general agreement between the members of this set of results, 
there are interesting and pertinent differences pertaining to spatial resolution. First, the 
secondary peak is not evident over the central Pacific Ocean at the 20°x30° grid scale, but 
appears over this region at the 10°xl0° and 5°x5° grid scales. This same effect is found over the 
Indian ocean and the northwest Pacific ocean subtropical regions. In addition, there are 
differences associated with the different horizontal scales on the occurrence of secondary 
maxima in the central-west Indian ocean, the central-east Pacific ocean, and north Africa. 

In summarizing the effects of spatial resolution on the areal cover of the secondary mode 
(see Table 1), it is found that the percentage of low resolution grid positions (upper panel in Fig. 
5) indicating the presence of a secondary peak is 78% (i.e., 14 of 18 grids), whereas for the 
medium resolution grid array (middle panel) the percentage falls to 69% (25 of 36). and then 
down to 65% (188 of 288) for the high resolution grid array. Based on the middle panel PR 
result of Fig. 3 at the highest spatial resolution used in the study (i.e., 5°x5°). the percentage rises 
to 68% (783 of 1 152). This emphasizes the sensitivity to grid resolution in evaluating the areal 
coverage of the secondary peak in the framework of 8-year means, with the end result being a 
difference of some 13% between the smallest and largest areal cover values. 

Notably, when these areal cover sensitivities to spatial resolution are decomposed into 
differences for oceanic and continental regimes, the percentage differences become larger, which 
as Table 1 indicates are 30% for ocean and 16% for continents. The much larger difference for 
ocean is primarily due to the mixing of continental and oceanic secondary peaks at the lowest 
resolution. Therefore, whereas it is evident that the secondary peak in precipitation’s diurnal 
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cycle is a widespread phenomenon over the global tropics and sub-tropics, caution must be 
exercised when evaluating the quantitative cover factors for oceans and continents because of the 
sensitivities involved. [Similar analyses using the TMI and combined PR-TM1 precipitation 
datasets lead to almost identical percentage differences, emphasizing robustness of the results.] 

Based on analyzing the 20th-century climate simulations by the newest generation of 1 8 
coupled climate system models. Dai (2006) indicates that most climate models can reproduce 
general precipitation patterns and their basic annual variability, although some of these models 
continue to generate unrealistic double Inter-tropical Conversion Zone (ITCZ) patterns over the 
tropical eastern Pacific ocean — a persisting problem in climate modeling (Mechoso et al. 1995). 
In addition. Dai (2006) found that most models produce too much convective precipitation and 
too little stratiform precipitation, resulting from unrealistically strong coupling of tropical 
convection to local sea surface temperatures. Most relevant to the concerns in this study, 
analyses of the precipitation diurnal cycle suggests that most models start to precipitate too early 
and too frequently at reduced intensity, showing a single late morning - early afternoon peak 
over continents and a single LE-EM peak over oceans. Many of these features are different from 
the TRMM precipitation diurnal analysis presented in this study, as well as in the Yang and 
Smith (2006. 2007) and Yang et al. (2006) studies. 

It should come as no surprise that climate models are unable to realistically simulate the 
diurnal variations of precipitation. Their spatial resolutions are too low vis-a-vis precipitation 
physics to properly simulate the narrowly confined mass and water vapor convergence processes 
leading to convection or to meaningfully resolve the vertical overturnings of cloud systems that 
actually produce convective and stratiform precipitation fallout. Moreover, they do not have any 
type of meaningful microphysical parameterizations insofar as the initiation, growth, loss, phase 
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change, and vertical motion of either non-precipitating or precipitating hydrometeors. Thus, it is 
appropriate that the type of observational ly-based study such as presented here, concerning 
highly detailed properties of precipitation processes, should be used to guide and calibrate future 
improvements in the embedded water cycle formulations of AOGCMs. 

4.2 Convective and Stratiform Partitioning of Secondary' Mode 

We have demonstrated that the horizontal distribution of the primary diurnal mode of 
precipitation over both oceans and continents is little affected by spatial resolution. However, 
the spatial averaging scale imparts a significant impact on the horizontal distribution of the 
secondary mode. Yang and Smith (2007) show that in an 8-year averaged global framework, the 
secondary mode is significant over continents but weak over oceans. They also find that when 
considering regional-seasonal scales the secondary mode becomes significant over oceans and 
even more so over continents, and the different modes are largely modulated by stratiform 
precipitation diurnal variability. In order to explain the possible cause of this secondary mode 
behavior, we analyze the diurnal cycles of precipitation’s convective and stratiform components. 

Figure 6 illustrates the horizontal distributions of the primary' and secondary diurnal 
cycles of precipitation from mean 8-year PR retrievals at the 20°x30° grid scale. By inspection it 
is evident that both convective and total precipitation exhibit the same oceanic - continental 
contrast in the phase of the primary diurnal peak. Also, the secondary mode of the convective 
component nearly mimics the secondary mode of total precipitation. Detailed inspection reveals 
that the phases of the dominant mode of both convective and total precipitation are 
approximately equivalent, however, there is a shift in phases between the secondary modes of 
convective and total precipitation. This shift is due to the contribution of the stratiform 
precipitation diurnal cycle. As seen in the lower panel of Fig. 6. an important feature is that the 
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dominant peak of stratiform precipitation's diurnal cycle over central Africa. Asia, Australia, and 
central South America occurs during the LE-EM period, while its secondary peak occurs during 
the MLA period at the same time as that of convective and total precipitation. 

Therefore, over continents a primary I. E-EM diurnal peak of stratiform precipitation is in 
phase with the secondary l.E-EM peak of total precipitation, while convective and total 
precipitation exhibit similar primary MLA diurnal peaks. T hese results indicate that afternoon 
convection forced by surface solar radiative heating reaching its maximum diurnal intensity is 
responsible for the dominant afternoon peak of the continental precipitation diurnal cycle. Le- 
the SRH mechanism. The development time of stratiform precipitation suggests that any one of 
a number of possible mechanisms such as the MTFPS, continental NRC-ERH. or perhaps even a 
SRC. DRC, or I..S-VM mechanism operating within a continental environment, might explain the 
LE-EM maximum in stratiform precipitation, which contributes to the secondary l.E-EM peak of 
total precipitation. In addition, the DRC mechanism operating within a continental environment 
might elevate nighttime precipitation while suppressing daytime precipitation — which is how the 
DRC mechanism operates. These processes of convective and stratiform systems draw attention 
to how multiple mechanisms might compete with one another for dominance, resulting in 
counterpoised LE-EM and MLA modes in the continental precipitation diurnal cycle. 

Over low latitude oceans, stratiform precipitation inevitably accompanies convective 
precipitation, especially for strong convective systems (Houze 1997). In fact, the contribution of 
stratiform precipitation to total accumulation is about the same as that of convective precipitation 
(Yang and Smith 2007). Fig. 6 actually shows that the dominant LE-EM peaks of convective 
and stratiform precipitation arc in excellent agreement, suggesting that whatever mechanism is 
producing the LE-EM mode, e.g.. one or more of the SRC. DRC. NRC-ERH, or LS-VM 
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mechanisms, is operating in synchronization with both convective and stratiform components. In 
addition, the secondary MLA mode, i.e.. due to the OSH mechanism, is effective for both 
convective and stratiform precipitation diurnal variability, although close inspection of Fig. 6 
suggests that the modulation of the ocean’s secondary mode amplitude stems mostly from 
fluctuations in the amplitude of the stratiform component. Similar analyses conducted at higher 
spatial resolutions lead to similar findings for the primary peak of the diurnal cycle, while the 
secondary peak exhibits somewhat more complex behavior (figure omitted). The secondary 
mode intermittently spreads out over the entire tropics and subtropics at the 5°x5° grid scale. 
Over oceans, this behavior suggests that the secondary MLA mode is more important at smaller 
horizontal scales. 

5. Harmonic Analysis of Diurnal Cycle 

5. / Diurnal Characteristics Revealed from Harmonic Analysis 

A Fourier harmonic decomposition scheme is applied to examine the diurnal variability 
of precipitation in order to elucidate the foremost amplitude and phase features while filtering 
high frequency noise. Because this approach suppresses spurious variations in a diurnal time 
series that might be detected as extrema in the explicit method, it is worthwhile evaluating its 
worth. Figure 7 compares the horizontal distribution of the precipitation diurnal cycle based on 
the explicit scheme applied to 5°x5° gridded PR-only retrievals over the 8-year ( 1998-2005) time 
period, to a counterpart result based on application of the Fourier decomposition scheme. The 
clock-face plots provide all essential details concerning amplitudes and phases, for not only the 
primary and secondary modes, but also for tertiary and quaternary modes. Colors and dock- 
hand pointing denote phase information while clock-hand length denotes amplitude information. 
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A strong contrast between primary harmonic diurnal phases over oceans and continents is the 
eye-catching feature in the upper panel of Fig. 7 (explicit scheme) as it is in the middle panel 
(Fourier decomposition scheme). The clear-cut contrast between the oceanic LE-HM phase and 
continental-like MLA phase along coastlines is actually sharper in the middle panel than in the 
upper panel. In addition, there are fewer spurious regions in the middle panel where the primary 
precipitation peaks are misidentified. A four-color clock-face scheme showing a more detailed 
representation of the phase interval of the primary mode is shown in the lower panel. Over 
oceans, frequent occurrence of late evening interval (0-6 MST) peaks is the dominant diurnal 
feature, with additional occurrences of morning interval (6-12 MST) peaks largely over coastal 
regions. Over continents, the occurrence of early evening interval (18-24 MST) peaks is the 
dominant feature over Africa and Australia, while the afternoon interval (12-18 MS) peaks are 
dominant over South America, North America and Asia. These more detailed diurnal features 
suggest that either different mechanisms are at work or singular mechanisms dominate but are 
dispersive in their timing of the maxima within different environments. Similar analyses at 
lower spatial resolution confirms that the underlying spatial distribution patterns of primary and 
secondary modes, as seen in Fig. 7, are not greatly affected by the use of a high resolution scale. 

The higher order harmonic modes conform to semi-diurnal, tri-diurnal. and quarlic- 
diurnal variations, and require an interpretation of where to position the actual peaks in terms of 
their phase angles in time. Figure 8 presents a comparison of the first three principal diurnal 
modes based on use of the explicit scheme in the upper panel, and the first two principal Fourier 
harmonics in the lower panel cast onto 20°x30° grids. The clock-face color exhibits a pre-noon 
half-day phase interval (blue) and post-noon half-day phase interval (green) of the primary 
diurnal mode. The secondary modes are shown by use of white inner clock-faces. It is obvious 
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that the horizontal distribution of the primary diurnal mode based on use of the explicit scheme is 
consistent with the horizontal distribution of the same mode based on use of the Fourier 
decomposition scheme, although close inspection of the two diagrams shows that small phase 
differences occur. More relevant are the greater and more frequent differences concerning both 
amplitudes and phases of the secondary mode between the two maxima identification schemes. 

The secondary mode will appear almost everywhere with Fourier analysis, simply 
because even if there is no secondary peak in actuality, almost invariably some power w ill leak 
into wave number 2. Even though it docs not satisfy the definition of the secondary mode over 
the northern part of South America and central South Africa, the secondary mode from the 
Fourier decomposition scheme is still apparent over these areas in comparison with the explicit 
scheme. For those grid positions where only a primary mode is found from the explicit scheme 
(such as over the central Pacific Ocean and the west Indian Ocean), a secondary mode is often 
identified from the Fourier analysis. The sources of the ambiguities associated with Fourier 
harmonic analysis of precipitation’s diurnal cycle are discussed in the next section. 

5.2 Ambiguities Associated with Harmonic Analysis 

The results presented in the previous section demonstrate that the Fourier decomposition 
approach is robust for the primary diurnal mode of precipitation, especially the phase of the 
dominant diurnal peak — as noted in. for example, the Collier and Bowman (2004) and Yang and 
Smith (2006) studies. However, the nature of Fourier harmonic analysis can lead to 
misinterpretation of both amplitudes and phases of different modes as they mix together, even 
the mixing of perfect sinusoidal modes. This can be especially problematic in seeking to identify 
and describe a secondary diurnal mode. Thus, ambiguities are inevitable when using Fourier 
analysis in analyzing precipitation’s diurnal variability. These ambiguities arise from a number 


17 



of sources. Of necessity, diurnal precipitation modes are: (a) generally truncated semi-diurnal 
modes, (b) represented as repeating cycles over the 24-hour diurnal period for wave number 1 
and beyond, and (c) not purely sinusoidal and thus dispersive in their power spectrum over a 
spread of frequencies when considering realistic diurnal modes in actual observations. 

Figure 9 demonstrates a case of Fourier decomposition for 2 purely sinusoidal, but 
truncated, diurnal modes with equal amplitude but shifted phase. The 2 peaks occur in early 
morning and late afternoon. The original components and combined time series for the double- 
peak cycle are illustrated in the left 3 panels. The right 4 panels present their successive 
reconstructed diurnal cycles using the mean and first 3 principal Fourier harmonics. Note that 
the primary mode is the 2 nd harmonic (wave number 2) which produces 2 diurnal peaks. Along 
with the mean value (wave number 0), the primary mode reproduces the basic diurnal properties 
but with a phase shift. Inclusion of the secondary mode (i.e.. 1 st harmonic, wave number l) 
produces a more realistic time series in which both amplitudes and phases are near their original 
values. However, it is evident that at least 3 harmonics are needed to reconstruct a sound 
representation of the original combined time series — based on the mixing of 2 perfect, but 
truncated, sinusoids. 

A similar analysis is conducted for another case in which the amplitude of the late 
afternoon peak is much smaller than that of the early morning peak (Figure 10). The primary 
mode (l sl harmonic, wave number 1) depicts nicely the phase of the principal diurnal peak, 
however, its amplitude is too small compared with the original value. Furthermore, the single 
primary harmonic cannot reproduce the secondary late afternoon peak. It is evident that the 
reconstructed diurnal cycle produced by the inclusion of the 2 nd harmonic (wave number 2) 
produces a realistic primary early morning peak in both amplitude and phase, as well as a 
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realistic secondary late afternoon peak in amplitude but with a phase shift. Thus, a 3 rd harmonic 
is needed to bring about a reconstructed diurnal cycle that is similar to the original. 

Figure 1 1 presents a final case, but now for a realistic observational situation, with a 
dominant peak at noon and a weaker secondary peak in late afternoon. Additional higher 
frequency fluctuations are also evident — which can be considered either real data properties, or 
more often noise due to retrieval error and/or under-sampling. Clearly, the harmonic approach 
can be used to eliminate the high frequency waves. In reconstruction, the primary mode (T 1 
harmonic, wave number 1 ) clearly reproduces the phase of the dominant peak, although the weak 
secondary peak is missed. Due to the small phase separation between the dominant primary peak 
and weak secondary peak, the 3 rd harmonic (wave number 3) represents the secondary mode. 
Using only the first two principal harmonics, the reconstructed diurnal cycle exhibits the basic 
features of the original cycle, however, with a discernible phase shift in the secondary mode. 

There are two main points to this analysis. First, when using the Fourier decomposition 
scheme, the dominant harmonic (along with the mean value), are able to reproduce the main 
features of the primary diurnal peak, particularly its phase. Second, the mean and first two 
principal harmonics generally can reproduce some, but not all. of the key features of the 
secondary peak. For some cases, the dominant harmonic cannot adequately reproduce the 
amplitude and/or phase features of the primary mode -- but for the TRMM data we have 
analyzed, these occurrences only arise ~6% of the time. Thus, since the diurnal cycle of 
precipitation often has a dominant primary peak and a weaker secondary peak, the principal 
Fourier harmonic is used to represent the main features of the strong primary mode, while the 
second harmonic is used to represent the main features of the weak secondary mode. In doing 
so. it must be recognized that a phase shift is often associated with the representation of the 
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secondary - mode, as well as possible amplitude errors for both modes. Note, the second 
harmonic is often needed to contribute to the amplitude behavior of the primary mode. Finally 
note, the harmonic approach can bring about a false secondary mode when only a single primary 
mode is present, in considering the case when a second harmonic is used to compensate for 
amplitude misrepresentation in the primary mode. 

Therefore, caution is always required in identifying secondary modes when using the 
Fourier decomposition scheme. Background knowledge of regional precipitation properties is 
always helpful in mitigating against ambiguities caused by applying harmonic analysis. In terms 
of our secondary mode analyses with TRMM precipitation data, the incidences of mis- 
interpretation associated with the Fourier decomposition scheme are relatively frequent, 
particularly for the higher spatial resolutions, contaminating the results some 45% of the time. 

6. Discussion and Conclusions 

This study has examined the diurnal variability of precipitation at various spatial and 
temporal scales based on 8-vear TRMM 1 MI. PR. and combined PR-TMI precipitation datasets. 
We find that these three distinct rain rate products produce consistent behavior insofar as the 
diurnal cycle properties of precipitation. The 8-year averaged diurnal cycles of oceanic and 
continental precipitation consistently exhibit primary maxima during the 3-6 MST and 15-18 
MST periods, respectively. Moreover, counterpart secondary maxima are strongly evident 
throughout the global tropics and sub-tropics, with seasonal variation being a very prominent 
feature of the secondary mode. These outstanding features of the diurnal behavior of oceanic 
and continental precipitation on a global scale are presented by use of a specially-designed 
graphic made up of gridded clock-faces with clock-hands, and coloring techniques to 
simultaneously display the various key parameters of harmonically complex precipitation 
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variation at the diurnal time scale. The secondary precipitation maxima generally consist of 
afternoon peaks over oceans and morning peaks over continents. Notably, this is the first study 
in which the widespread secondary mode has been identified and given special attention. 

it is important to recognize that based on our past and current TRMM data analyses, vve 
find that the secondary maxima are largely produced and modulated by stratiform precipitation 
(see Yang and Smith 2007). Although we have shown that the underlying properties of the 
precipitation diurnal cycle are influenced by the spatial and temporal resolutions at which the 
data are analyzed, we find that there is a great deal of consistency in how oceanic and continental 
regions contrast with one another, regardless of the underlying resolutions. For example, the 
impact of spatial resolution is generally not significant concerning the primary diurnal mode, 
although it can be significant in identifying the areal coverage of the secondary diurnal mode. 
Most importantly, the global distribution of the precipitation diurnal cycle exhibits considerable 
seasonal variations — largely due to the occurrence of secondary diurnal modes. 

The Fourier decomposition scheme does a credible job in identifying the main amplitude 
and phase features of the primary mode with the advantage of eliminating spurious high 
frequency fluctuations in the original diurnal time series. We find that this approach is at its best 
when representing the phase of the primary diurnal mode. Nonetheless, at times, secondary, 
tertiary, and even quaternary modes are needed in reproducing realistic amplitudes of the 
primary mode. Thus, the Fourier decomposition scheme can produce isolated ambiguities vis-a- 
vis the interpretation of the primary diurnal mode, depending on how the scheme is applied. 

Furthermore, we have demonstrated how and why the Fourier decomposition scheme 
creatse ambiguities in representing the amplitude and phase properties of the secondary diurnal 
mode because of the up-frequency mixing of power from the first principal mode into the higher 
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order modes. We have shown that secondary diurnal modes over some regions derived from the 
Fourier scheme applied to TRMM data — are simply not real. Therefore, caution is in order 
when interpreting the amplitudes and phases of the secondary, tertiary, and quaternary modes, 
apart from any coherency at the global scale, when using Fourier analysis — noting that 
ambiguities in representing the secondary mode are unavoidable. Nonetheless, surface 
precipitation normally exhibits fairly distinct primary and secondary diurnal peaks such that the 
Fourier decomposition scheme is often applicable. In any case, background knowledge of the 
pertinent meteorological variables controlling precipitation, as well as understanding how the 
precipitation variables behave in response to the relevant diurnal mechanisms at work, are 
always helpful in analyzing for secondary diurnal modes when using Fourier decomposition. 
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Table 1: Areal coverage of secondary peak in diurnal precipitation cycle (given in %) at various 
spatial resolutions. [Results based on TRMM PR measurements.] 



5°x5° 

10°xl0° 

20°x30° 

20°x60° 

Continental 

Regime 

68 

64 

85 

80 

Oceanic 

Regime 

69 

65 

55 

85 

Total 

68 

65 

69 

78 


29 









List of Figures 

Figure t: Mean diurnal cycles of oceanic (left panel) and continental (right panel) precipitation 
from eight years (1998-2005) of TRMM satellite data based on TMI (2al2). PR (2a25). and 
combined PR-TMI (2b31) TRMM standard algorithm retrievals. [Key to both diagrams 
given in left panel.] 

Figure 2: Mean diurnal cycles of seasonal precipitation from PR algorithm retrievals for eight 
years (1998-2005) over oceans (left panels) and continents (right panels). Total, convective, 
and stratiform precipitation classifications are denoted by T, C. and S. respectively. [Key to 
all diagrams given in upper-left panel.] 

Figure 3: Horizontal distributions of diurnal precipitation cycles based on TRMM satellite data 
averaged over 8-year period (1998-2005) and analyzed on 5°x5° spatial grids — using clock- 
face diagrams. Upper, middle, and lower panels provide results for TMI, PR. and combined 
PR-TMI algorithms, respectively. North, east, south, and west points of given clock-face 
refer to 00. 06. 12, and 18 MST phases, respectively. Blue clock-face indicates that primary 
mode occurs some time during pre-noon half-day (0-12 MST) phase interval, while green 
clock-face indicates that primary mode occurs some time during post-noon half-day (12-24 
MST) phase interval. Black upper half-circle of each clock-face denotes nighttime period, 
while white lower half-circle denotes daytime period. White inner clock-face highlights grid 
position where secondary diurnal maximum is clearly evident. Red clock-hand points to 
numeric phase of primary mode, green hand of secondary mode, and blue hand of tertiary 
mode — over 0-24 MST interval. Amplitude of individual mode is denoted by length of 
associated clock-hand, in which four embedded rings indicate magnitudes (on log scale) of 

0.15. 0.3. 1.5. and 3 mm day '. 
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Figure 4: Same as in Fig. 3. except that results are analyzed on 20°x30° spatial grids such that 
rows represent 10-30°S. 10°S-10°N, and 10-30°N latitude belts. 

Figure 5: Same as in Fig. 3, except for only PR algorithm with separate panels illustrating 
results on 20°x60° (upper). 20°x30° (middle), and 10°xl0° (lower) grids. 

Figure 6: Same as in Fig. 4. except for only PR algorithm with separate panels illustrating 
results for total (upper), convective (middle), and stratiform (lower) precipitation. 

Figure 7: Comparison of horizontal distributions of PR algorithm- based precipitation diurnal 
cycles involving two different maxima identification schemes applied in 8-year averaging 
framework (1998-2005) on 5°x5° grids — using clock-face diagrams. Upper panel shows 
dominant maxima based on "explicit" scheme. Middle and lower panels show first four 
dominant maxima based on “Fourier decomposition" scheme — using two different clock- 
face color formats for indicating modal phase intervals ( l sl of these formats is same as used 
for upper panel). Two-color clock-face format is used for upper and middle panels, in which 
blue clock-face denotes that primary mode occurs some time during pre-noon half-day (0-12 
MST) phase interval, while green clock-face denotes that primary mode occurs some time 
during post-noon half-day (12-24 MST) phase interval. Four-color clock-face format is used 
for lower panel, in which blue, light blue, green, and light green faces denote that primary 
mode occurs some time during late evening (0-6 MST). morning (6-12 MST). afternoon ( 12- 
18 MST), or early evening (18-24 MST) phase interval, respectively. Red clock-hand 
indicates numeric phase of primary harmonic, green hand of secondary harmonic, blue hand 
of tertiary harmonic, and black hand of quaternary harmonic — over 0-24 MST interval. 
Amplitude of individual mode is denoted bv length of associated clock-hand, in w hich four 
embedded rings indicate magnitudes (on log scale) of 0.15, 0.3, 1.5, and 3 mm day . 
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Figure 8: Same as in Fig. 7. except for: (a) application on only two 20°x30° grids, (b) use of 
only two-color clock-face format, (c) extraction of only first three principal modes for 
explicit scheme, (d) extraction of only first two principal harmonics for Fourier 
decomposition scheme, and (e) use of white inner clock-face to denote grid position where 
secondary diurnal maximum is clearly evident. 

Figure 9: Fourier decomposition of two purely sinusoidal, but truncated, equal-amplitude and 
phase-shifted diurnal modes combined into two-peak diurnal cycle. Three left-hand panels 
show two separate diurnal modes along with combined diurnal cycle. Four right-hand panels 
show 7 successive reconstructions of combined diurnal cycle using mean value (wave number 
0) and first three principal Fourier harmonics (wave numbers 1 - 3). 

Figure 10: Same as in Fig. 9. except that two diurnal modes are unequal-amplitude. 

Figure 11: Same as in Fig. 9. except that two diurnal modes are realistic (non-sinusoidal and 
unequal-amplitude). 
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Figure 1: Mean diurnal cycles of oceanic (left panel) and continental (right panel) precipitation 
from eight years (1998-2005) of TRMM satellite data based on TMI (2al2). PR (2a25). and 
combined PR-TMI (2b31) TRMM standard algorithm retrievals. [Key to both diagrams 
given in left panel.] 
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Figure 2: Mean diurnal cycles of seasonal precipitation from PR algorithm retrievals for eight 
years (1998-2005) over oceans (left panels) and continents (right panels). Total, convective, 
and stratiform precipitation classifications are denoted by T. C. and S. respectively. [Key to 
all diagrams given in upper-left panel.) 
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Figure 3: Horizontal distributions of diurnal precipitation cycles based on TRMM satellite data 
averaged over 8-year period (1998-2005) and analyzed on 5°x5° spatial grids — using clock- 
face diagrams. Upper, middle, and lower panels provide results for TMI, PR, and combined 
PR-TM1 algorithms, respectively. North, east, south, and west points of given clock-face 
refer to 00. 06. 12. and 18 MST phases, respectively. Blue clock-face indicates that primary 
mode occurs some time during pre-noon half-day (0-12 MST) phase interval, while green 
clock-face indicates that primary mode occurs some time during post-noon half-day (12-24 
MST) phase interval. Black upper half-circle of each clock-face denotes nighttime period, 
while white lower half-circle denotes daytime period. White inner clock-face highlights grid 
position where secondary diurnal maximum is clearly evident. Red clock-hand points to 
numeric phase of primary mode, green hand of secondary mode, and blue hand of tertiary 
mode — over 0-24 MST interval. Amplitude of individual mode is denoted by length of 
associated clock-hand, in which four embedded rings indicate magnitudes (on log scale) of 
0.15, 0.3, 1.5, and 3 mm day '. 
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Figure 4: Same as in Fig. 3. except that results are analyzed on 20°x30° spatial grids such that 
rows represent 10-30°S, 10°S-10°N. and 10-30°N latitude belts. 
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Figure 5: Same as in Fig. 3. except for only PR algorithm with separate panels illustrating 
results on 20°x60° (upper), 20°x30° (middle), and 10°xl0° (lower) grids. 
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Figure 6: Same as in Fig. 4. except for only PR algorithm with separate panels illustrating 
results for total (upper), convective (middle), and stratiform (lower) precipitation. 
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Figure 7: Comparison of horizontal distributions of PR algorithm-based precipitation diurnal 
cycles involving two different maxima identification schemes applied in 8-year averaging 
framework (1998-2005) on 5°x5° grids - using clock-face diagrams. Upper panel shows 
dominant maxima based on '‘explicit" scheme. Middle and lower panels show first four 
dominant maxima based on "Fourier decomposition” scheme — using two different clock- 
face color formats for indicating modal phase intervals (l sl of these formats is same as used 
for upper panel). Two-color clock-face format is used for upper and middle panels, in which 
blue clock-face denotes that primary mode occurs some time during pre-noon half-day (0-12 
MST) phase interval, while green clock-face denotes that primary mode occurs some time 
during post-noon half-day (12-24 MST) phase interval. Four-color clock-face format is used 
for lower panel, in which blue, light blue, green, and light green faces denote that primary 
mode occurs some time during late evening (0-6 MST), morning (6-12 MST), afternoon (12- 
18 MST), or early evening (18-24 MST) phase interval, respectively. Red clock-hand 
indicates numeric phase of primary harmonic, green hand of secondary harmonic, blue hand 
of tertiary harmonic, and black hand of quaternary harmonic — over 0-24 MST interval. 
Amplitude of individual mode is denoted by length of associated clock-hand, in which four 
embedded rings indicate magnitudes (on log scale) of 0. 1 5. 0.3. 1 .5, and 3 mm day . 
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Figure 8: Same as in Fig. 7. except For: (a) application on only two 20°x3()° grids, (b) use of 
only two-color clock-face format, (c) extraction of only first three principal modes for 
explicit scheme, (d) extraction of only first two principal harmonics for Fourier 
decomposition scheme, and (e) use of white inner clock-face to denote grid position where 
secondary' diurnal maximum is clearly evident. 


40 


Explicit Fourier 


I Diurnal Component 1 1 


/ 

J 


/ osoolst 


\ 


Local SoUsr Time 


I Diurnal Component 2 { 


! I 


-'l 3 
I 
! 

| 2 
5 

\ I 

I 

j 

5 « 

5 

f< -i 
3 
2 
1 


Local Solar Time 


Actual Diurnal Cycle 



0 6 12 18 24 

Local Solar Time 


3 


2 


I 


0 


1 


Individual Harmonic of Oiumal Cycle (ordered by contribution to total variance) 
Progressively Reconstructed Diurnal Cycle 
Actual (original) Diurnal Cycle 


Q. 

E 

< 


Fourier Component 1 (wavenumber 0) 


Mean Value 


6 12 18 
Local Solar Tima 



6 12 18 
Local Solar Time 


2 

Q- 

E 

< 


fj- -------- 3 

Fourier Component 2 (wavenumber 2) 

Fourier Component 4 (wavenumber 3) 

i 2 
. 

/ • \ / \j 

\y . / \ / \ m 

\ : o 

- 

A A 1 
/ ^ \ 


o 6 12 18 24 0 6 12 18 24 

Local Solar Time Local Solar Time 


Figure 9: Fourier decomposition of two purely sinusoidal, but truncated, equal -amplitude and 
phase-shifted diurnal modes combined into two-peak diurnal cycle. Three left-hand panels 
show two separate diurnal modes along with combined diurnal cycle. Four right-hand panels 
show successive reconstructions of combined diurnal cycle using mean value (wave number 
0) and first three principal Fourier harmonics (wave numbers 1 - 3). 


41 








opnmduiv apnitiduiv 


Individual Harmonic of Oiurnal Cycle (ordered by contribution to total variance } 
Progressively Reconstructed Diurnal Cycle 
Actual (original) Diurnal Cycle 


| Diurnal Component 1 
0800 LST 

j \ 


* U |* ~ M 

local Solar Time 


| Diurnal ComporvenTT 


1700 LST 



» I II » 14 

Local Solar Tim© 


Actual Diurnal Cycle 


6 12 18 24 

Local Solar Time 


© 

~c 

3 

a 

I 


Fourier Component 1 (wavenumber 0) 

Fourier Component 3 (wavenumber 2} 

i > 

i ; 
1 1 | 

i. i 
: | 

\ | 

1 Y 

» / 

\ { 

A 

l 


» 6 12 IS 24 0 6 12 IS 24 

Local Solar Time Local Solar Time 



1 3 : 


Fourier Component 2 (wavenumber 1) 


Fourier Component 4 (wavenumber 3) 


2 



0 6 12 18 24 1 II 6 12 IS 24 

Local Solar Time Local Solar Time 


Figure 10: Same as in Fig. 9. except that two diurnal modes are unequal-amplitude 
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Figure 11: Same as in Fig. 9. except that two diurnal modes are realistic (non-sinusoidal and 
unequal-amplitude). 
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